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Abstract 

An energy-minimal simulation is proposed to study the patterns and mechanical properties 
of elastically crumpled wires in two dimensions. We varied the bending rigidity and stretching 
modulus to measure the energy allocation, size-mass exponent, and the stiffness exponent. The 
mass exponent is shown to be universal at value Dm = 1-33. We also found that the stiffness 
exponent a = —0.25 is universal, but varies with the plasticity parameters s and 9p. These 
numerical findings agree excellently with the experimental results. 

PACS numbers: 62.20.-x, 05.90.+m, 89.75.Da 
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I. INTRODUCTION 



The process of crumpling is everywhere in nature and human activities, including the 
formation of mountains and valleys in tectonics[l|, packing of DNA strands in viruses^, 
car wreckage after an accident sl , or the noisy food wraps that drive us nut in the theater, 
etc.. In spite of its ubiquity, the mechanism behind many of its properties has remained 
unclear ^, |5| . For instance, how does the labyrinthian internal structure evolve such that it 
can withstand extraordinary pressure while more than 80% of its interior remains vacant? 
Also, why is it that there exists a power law between the external force and the sphere radius 
with an exponent that varies with material 6| but, otherwise, is insensitive to the thickness 
and size of the thin sheet? 



Scientists have studied crumpled wires (CW) both theoretically 



9| and 



experimentally 1^, [ll|, [l^ . As is shown in FiglH the wire is smooth and rid of the com- 
plicated ridges and vertices, which is different from a crumpled sheet [isl. In the last few 
years, more and more interesting properties of CW have emerged. For instance, Donato et 
al. found a scaling law in the size-mass ..ationQ, while Offered the a.„c,a.ed 

morphological phase diagram and reported a power law for the number of loops with dif- 
ferent exponents in each morphology 7| • However, there are still few quantitative studies of 
CW either on the energy aspect or macroscopic properties. Furthermore, the former authors 
all worked on stuffing the wire into a fixed two-dimensional cavity. It is not clear whether 
this shares the same properties as by increasing the strength of the confining potential while 
fixing the wire length. Aside from this, another motivation for us is to check how sensitive is 
the CW on the specific form and relative amplitude of stretching and bending energies. In 



this report, simulations are carried out to study the mechanica 
semifiexible polymer model 9|] and the minimal energy model 



properties of CW under the 



14l |. To prevent the entropy 



from dominating the statistical behavior [15[, the wire is chosen to be of medium length. 



II. MODEL AND COMPUTATIONAL METHOD 



In our simulations, the optimal shape of a wire is obtained by minimizing the total 
energy of CW. We begin with a wire in a weak external field, and gradually increase the 
field strength. At each stage with a certain field strength, we recursively minimize the 
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FIG. 1: A crumpled copper wire in a two-dimentional cavity. This photo is reproduced by courtesy 
of Donato et aL|10|. 

energy and move to the next stage with a larger external field after an equilibrium state has 
been achieved. This process is carried out throughout the whole simulation until the wire 
structure no longer changes with the increasing external field. The wire consists of A'^ number 
of monomers and is of length Lq = 1. The total energy of a CW includes bending energy 
-Efe, stretching energy E^, and external potential -Eext- In addition, a hard-core potential E^f. 
is incorporated to keep the wire from self-crossing. Explicitly, the energy function can be 
written as: 

N-l , ^ t. 

Etot = E + E [|(A/.)^ + Xr^] + E,, (1) 

1=1 i=l 

where 6i is the bending angle between two successive segments and = YI /{Lq/N) is the 
bending rigidity which is proportional to the Young's modulus Y and second moment of 
inertia I. The stretching modulus kg = Y A/ {Lq/N) depends on the wire cross section A, 
and AZj is the deviation of the z-th segment from its equilibrium length Lq/N . The distance 
Ti of the z-th monomer is measured from the the center of the external parabolic potential. 
Since most research interests fall into cases with insignificant extension, we mainly focused 
on conditions that ki, <^ k^. Typical parameters for simulations are listed in Tab. [H The 
steady state of CW is determined by minimizing its total energy with Powell's algorithm. 
To make sure that the local minimum we found was not a special case, we adopted a random 
set of searching directions in the Powell's algorithm and repeated the simulations for ten 
times to get the averaged results. Most of the time, CW are randomly folded rather than 
forming a perfect spiral - the apparent global minimum [l^. This indicates that those folded 
configurations are metastable states. The modified Powell's minimization is equivalent to 
the Monte Carlo simulation in this aspect. 
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TABLE I: The value of parameters used in simulations 



Parameter Range of values 

Bending rigidity kb 1 to 10 

Stretching modulus ks 10^ to 10^ 

External field A 1 to 2 x 10^ 

Number of segments 100 

Wire length Lq 1 

Width of hard core 10~^ 
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FIG. 2: A sequence of simulation results for fc^ = 5 and fc^ = 5 x 10^. The value of A for (a), (b), 
(c), and (d) are 70,200,300, and 1000 respectively. It can be seen that when we increase A, no 
sharp edge occurs and the shape of CW remains smooth. Those local structures start to contact 
each other when A > 1000. More turning points are formed if the external field increases. 



III. RESULTS AND DISCUSSION 



A. shape and size-mass exponent 

The loops of CW in FigsiU and [2] resemble the boundary of a water-drop- like structure. 
Compare with the simulations for a semi-flexible polymer [3| which explicitly put in finite 
temperatures and consequently resulting in many small fluctuations, the configuration of 
CW is characteristically smooth, and so we exclude the temperature effect in this simulation. 
We believe the smoothness is credited to the stretching force which causes the relaxation 
of roughness, analogous to the surface tension which disfavors a kinky surface. Because 
the stretching energy is found to be insignificant when kf, <^ k^, we can neglect it when 
determining the mathematical form of the configuration. For simplicity and without the 
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n k^=3kj=2.4xlo' 
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Radius of circle R 



FIG. 3: The mass-size relation of CW for A=1000. The radius of the circle is rescaled by Rg- 
Apparently, all mass will already be included in the circle when R= 1. The rescaling process only 
shifts the data on a Log-Log plot and does not affect the value of exponent Dm = 1-33. Since we 
set the origin of the circle at the center of mass of the CW, there is a bias in the data at small R. 
This bias is avoided by fitting only from i?=0.5 to R=l as shown by the solid straight line[lo|. 

loss of generality, we consider the case oi N ^ oo, i.e., the continuous limit. Eq.([T]) then 
can be re-written as: 



7-1* 

^tot 



2 ^ds' 



ds 



(2) 



The search for shape function r{9) with minimal energy is done by introducing the variational 
method. Fix both r(s = 0) and r(s = L) and this total energy is identical to the classical 
action once we identify the arc length parameter s as the time. The Euler-Lagrange equation 



of motion immediately gives out the shape function as r'^ \fkh[\ = sin 26', where the term 
is not hard to expect by the dimensional analysis in retrospect. This matches the pattern in 
the initial stage of the crumpling observed experimentally lOf] as well as in our simulations. 
However, when the loops start to touch each other as the crumpling proceeds further, the 
excluded volume or steric interaction 



15| is expected to squeeze and distort the shape of 



loops. 



n 



The mass function M{R) is a common property[10|] to characterize the configuration of 
CW. It is defined as the total mass encompassed within a circular area of radius R centering 
at the CW's center of mass: 



M{R) =Y,d{R- n)5{e{n - R), ^(r.+i - R)) 



i=l 



N 



(3) 



-6ie{n - R),eiR - r,+i))- fin, r,+0 
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where 6{x) denotes the step function, S{x,y) stands for the Kronecker delta-function which 
equals 1 when x = y and zero otherwise, and f{ri,ri^i) is the portion of length between 
and rj_|_i within the circle. After performing a series of simulations, the mass-size relation is 
computed and shown in Fig 13] for six distinct combinations of and kg. We can see that the 
mass of CW grows in a power-law fashion when the radius increases and, disregarding jl^ 
the fluctuations at small external field, the growing exponent saturates at Dm = 1.33 for A > 
100, see Figini The saturating behavior is in accordance with Donato et a/.'s conjecture flo| 
that in a loose-packing situation the corresponding mass-size exponent has insignificant 
dependence on how the wires are injected. And the value of Dm = 1-33 is also consistent 
with Donato et a/.'s result. 



B. Energy allocation 

We characterize the allocation of elastic energies in CW by the ratio '-f=Es/Ei). Data in 
Fig. m establish a scaling relation between 7 and the three parameters in the Hamiltonian 
as: 

7 ~ A^' kl^ kl^ (4) 

where ei = 0.55, = 2.6, and = —3.3. The numerical results can be deduced analytically. 
These exponents can be obtained by doing a dimensional analysis on the bending energy and 
the external potential in Eq.([2]), kb/i ~ Xi^, which gives the characteristic size for each loop 
as £ ~ (kb/Xy^^. Multiplying the bending energy stored in each loop, kb/i, by the number 
of loops L/i enable us to estimate the bending energy as kbL/i^. In a similar fashion, the 
stretching energy can be written as Eg ~ ks{L — Lq)"^. Minimizing the sum of these two 
energies with respect to L gives 

^ ^ VhX 

ks [Lq - y/hX/ks] 

The value of ei = 0.5 matches the simulation result. Although the exact magnitudes are 
not the same, the signs of and eg and their relative magnitude are still captured. One 
possible explanation for the discrepancy is our failure to include the hard-core potential in 



the analytic arguments, which has been known to be difficult ISj but crucial for a realistic 
polymer. 
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FIG. 4: The value of 7 is plotted against several parameters. For (a) and (b), 7 grows with the 
strength of external field in a power-law fashion with exponent e/ = 0.55. In (a) the stretching 
modulus is set at 2 x 10^, and in (b) the bending rigidity is 6. Panels (c) and (d) depict 7 versus 
ki, and kg respectively, and they share the same denotation of data. Notice that the power-law can 
only be found for strong external field (A > 500). The exponents of the power-law are determined 
to be eb = 2.6 and = —3.3 
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FIG. 5: The shrinkage of Rq versus the strength of confining potential for CW with different 
rigidity and modulus, kh and kg. All data collapse to the parameters kb = 8,ks = 2.4 x 10^ with 
the scaling Rg ~ {k^ksY'"^- The solid line is the fitting result Rg ^ A"*^'^^. 

C. Stiffness and plasticity 



How to estimate the resistance of a crumpled structure is a well-known difficult problem. 
For instance, what is the main reason why such a loose crumpled structure can be so hard? 
As a starting point, people has studied the power law between the structure size and the 
applted force both experimentally and tlteoret,callyB,Q . Our simulation result in FigjS] 



showed that the radius of gyration Rg 



of CW decreases with A in a power-law form. 
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To prevent the structure from becoming too complicated, we limit the field strength A at 
2000 to optimize the resolution and accuracy. The data are found to coUape into the scaling 
relation Rg ~ X°'{ki)ksY-^, which suggests that the stretching energy not only correlates those 
monomers, but also contributes to the stiffness. Intuitively, the dependence of kg is expected 
to drop out instead of diverging as it approaches infinity. The reason why we obtained such 
a scaling relation is a consequence of the discrete model since the length of monomers is 
bound to be compressed as they all scramble to be near the center of the external potential. 

In FigJSl a = —0.25 is determined to be independent of kb and kg and, consequently, the 
thickness and Young's modulus of the wire. This result is the same as the predictions made 

n 

for an elastic thin sheet [18|]. It can be derived analytically by minimizing with respect to R 
the sum of the external potential, XR^Lq, and the total repulsive interaction, I3{Lq/ R'^Y -R^ 
where /5 is a measure of the hard-core potential, L^/R^ is the monomer density, and the 
other R^ factor comes from integrating over the cavity. 

Since Stoop et al. has demonstrated the importance of plasticity for CW, we also include 
it in our simulation. The bending rigidity kb is revised with a linear approximation, skb{^ < 
s < 1) in the plastic regime, so that the bending energy becomes Eb = kb9l/2 + skb(6^ — 
9p) /2 + kb9p{l — s){9 — Op) when 6 is larger than the yield threshold angle 6p\i\. We analyzed 
the \-Rg relation of CW for s = 0.05,0.1,0.2,0.4, and 0.8, and it behaves very different 
from that of the elastic one. In FiglHl the value of a decays with the bending rigidity in 
the plastic regime which is consistent to the experimental findings in j^. It also decays 
exponentially with the yield threshold 6p. When large deformations cost less energy due to 
the plasticity, they appear in abundance and cause many vertex-like structures which makes 
the CW stiffer. 



IV. CONCLUSION 

In conclusion, we introduced a minimal one-dimensional model to simulate the crumpled 
wires. The simulation results agree excellently with the experimental observation, which 



10|. We 



include the pattern of loops, the universal value of mass-size exponent Dm = 1-33 
checked the exponent a that characterizes the power-law relation between the crumping force 
and the CW radius to be independent of the bending rigidity and the stretching modulus 



which is consist to 18|. Plasticity is found to suppress a which is in the right trend as the 
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FIG. 6: The X-Rg relation of the CWs{{kb, kg) = (5, 5 x 10^)) with plasticity for different s which 
are denoted in the figure. The value of the exponents a are -0.21, -0.24, and -0.29 for s=0.4, 0.2, 
and 0.05 respectively. The inset illustrates that a decays with 9p exponentially and eventually 
saturates at value of 0.26 for s = 0.4. 

recent experiment [gI on crumpled sheets. 
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